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I ^ ■ ABSTRACT 

. , We discuss the question to what degree the geometrical structure and the matter content of the 

' universe at scales exceeding the present Hubble horizon is constrained by cosmological observations, 

^_^ , in particular by measurements of the cosmic microwave background radiation. For an answer, a 

^ ' simple formalism is described, which goes back to a paper by Kristian and Sachs in 1966. 

^: 
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O ! 1. INTRODUCTION 

o\ . 

0^ ■ It is usually assumed that the observable universe can be described by one of the 

^^ I homogeneous-isotropic Friedman-Robertson- Walker models, at least up to scales of 

O ■■ the order of the horizon Xh ■ The structure of the universe at larger scales is not 

known. The inflationary scenario suggests, that the observable universe is part of 
an inflated region of size Xj , possibly much larger than Xh ■ Outside the inflated and 
"ti I hence nearly homogeneous region the spacetime geometry might be inhomogeneous 

C^ ■ and anisotropic. Could structures outside the horizon by means of their tidal fields 

^ I have some imprint on the observationally accessible part of the universe (shortly: the 

•'^ ■ local universe), in particular, may they influence the cosmic microwave background 

r> I radiation? 

S ■ Twenty years ago L. Grishchuk and Ya. B. Zeldovich (1978) (see also Turner 

■ ■ ■ ' 1991, Grishchuk 1992, Kashhnsky et al. 1994, Lyth 1995, Garcia-Belhdo et al. 

1995) have discussed a similar question. They found that the CMB perturbations 
6T/T resulting from the growing mode of density or gravitational wave perturba- 
tions on a scale A and corresponding to a metric perturbation h are of the order 
ST/T ^ ft,(A///A)^. Since one expects h ^ 1 for regions A > A/, the COBE obser- 
vations lead to A/ > 10'^ Xh, thus the "homogeneous region" must be much larger 
than the Hubble distance, independent of any inflationary hypothesis. This argu- 
ment does not exclude structures exceeding Xh only slightly, if they correspond to 
metric amplitudes h much smaller than 1. Also other deviations from FRW models 
may be allowed to a certain degree. For instance, if one restricts the geometry to 
homogeneous Bianchi models, one is able to obtain from the CMB measurements 
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upper limits for the vorticity and shear of these models (Hawking 1969, Collins and 
Hawking 1973, Bunn et al. 1996, Kogut et al. 1997). But constraining the geometry 
to a small class of simple models may be misleading, since the real universe might 
have a much more complicated geometrical structure and matter content at very 
large distances. In view of the remarkable improvements expected for future CMB 
observations with the MAP and PLANCK missions it is of interest to see how much 
of this structure could be detected or constrained. A first step is a mathematical 
framework dealing with sufficiently general cosmological models. We describe a 
very simple proposal, which tackles the problem, but needs further elaboration. 

3. KRISTIAN-SACHS APPROACH 

J. Kristian and R.K. Sachs (1966) made in their classical paper only few assump- 
tions: (i) The universe is described by a Riemannian spacetime with slowly vary- 
ing metric tensor; (ii) light travels along null geodesies and obeys the usual area- 
intensity law; (iii) the gravitational field is related to the matter by the Einstein 
field equation for dust; (iv) all quantities of interest can be expanded in a power 
series around here-and-now as origin. In their treatment it not necessary to adopt 
a specific cosmological model a priori, the results are complete and general (apart 
from the assumption of dust, which can quickly be generalized). Ignoring global 
properties (or rather assuming that the scale of a possible compactification of the 
universe is larger than all other scales) made the calculation manageable, and last 
not least all introduced quantities are measurable in principle. 

With a power series expansion around here-and-now one can hardly treat evo- 
lutionary problems in cosmology. In particular, the CMB fits not easily into the 
framework. We remove this shortcoming by treating space and time coordinates on 
a different footing. Temporal variations are seen to appear with much higher ampli- 
tude than spatial variations - at least if one averages over the inevitable small-scale 
fluctuations as origin of galactic structures. This leads to the idea to carry out the 
power-law expansion only within spatial hyper-surfaces orthogonally intersecting 
the observers world line. The Kristian-Sachs expansion coefficients then become 
time-dependent. Actually, we reduce the generality by adopting the (anthropo- 
morphic) assumption, that space-time reduces along the observers world line to the 
homogeneous-isotropic FRW models. In geometrical terms, the observers world line 
L is assumed as geodesic and shear-free, generally in contrast to the world-lines of 
neighbouring observers. The method remembers to some degree the covariant ap- 
proach put forward by G. Ellis and his coworkers in numerous publications, but 
contrary to their work we try to fix the coordinate system as much as possible in 
order to reduce the number of coefficients needed to describe deviations from FRW 
models. 

3. COORDINATE SYSTEM 

Basic to the method is the construction of a coordinate system, which is adapted 



to the matter distribution. We represent the cosmic matter as relativistic fluid 
flow with V^ as the four- velocity vector tangent to the flow hues {V^Vfj. = —1, 
our notation foUows the book by Misner, Thorne and Wheeler 1973). Comoving 
coordinates are introduced such that V^ = Sq. The flow lines are given by x* = 
const (i — 1, 2, 3). The time coordinate t is not completely fixed by its coincidence 
with the proper time on the flow lines, and also the spatial coordinates are subject 
to arbitrary time-independent changes, t = t + /(x*), x^ = a;*(x'^). 

Let the particular flow line x* = be the observers world line L. We demand 
that on and near L a FRW solution (with zero spatial curvature, see below) is a 
good approximation. It is suflicient to assume that L (and only this world line) is 
geodesic and shear- free. This is a weak form of reducing to FRW along L, since 
it is compatible with a non-vanishing Weyl tensor and a nonzero vorticity on L. 
We may then represent the metric tensor near L as a Taylor series in the spatial 
coordinates x': 

500 = -1 (1) 

5oj = kkx'' + kkix''x'' +kkinix''x''x"'' + ... (2) 

gik = a^Sik + hikix^ + htkimX^x"^ + hikimnx'-x"^x"' + ... (3) 

where a(t) is the Friedman scale factor, the other expansion coefficients are as well 
time functions in general. The coordinate transformations preserving this form of 
the metric may also be expanded with constant coefficients, which later serve to 
reduce initial values of the metric expansion coefficient: 

t = t + akix'^x'- +aki„ix''x'-x"' + ..., (4) 

x' = x' + bikix'^J + hiklmX^X^x"^ + .... (5) 

Only in a linear approximation the terms added to the FRW metric can uniquely 
be separated into scalar, vector and tensor (gravitational wave) perturbations, we 
shall not attempt this here. - We have assumed a flat FRW metric along L, since 
the spatial curvature of the FRW models is hidden in the expansion coefficients of 
Eqn (3). 

4. MATTER CONGRUENCE 

Using the decomposition of the velocity gradient 

V^^.u = ojt^u + cr^^ + 0(5^^ + V^,V^)/2, - K^M, (6) 

we can expand acceleration, shear, expansion and vorticity around L: 

A, = {kk+kkix^)x^ +o{i), (7) 

CTkl = -{hklr - -5klhssr)x'' +o{2), (8) 

e = 3^ + ^|(%^)x'+o(2), (9) 

a 2dt a'^ 



i^ki = -^{hi - hk) + ihir - kkr)x'' + o(2), (10) 

together with Aq = 0, ctq/^ = 0, LUok = 0. The lowest terms in this expansion 
transform under (3), (4) as 

hiit) = hi + ^aki, (11) 

hklm{t) — hklm{t) — 2a (bklm+blkm)- (12) 

Since aki = aj/t, the vorticity uj^i on L is not affected by coordinate transformations. 
Similar simple transformation rules hold for the higher-order coefficients. For the 
matter tensor we assume a perfect fluid 

7)iiy = Vfj,Vi,{fi + p) + pgpiu (13) 

with an equation of state p = w^ {w — const). The conservation equations T.^'^ = 
can partly be integrated without expanding /i and p in a Taylor series. One obtains 

^i = ™(x*)((l + 5ofc5oa'')rfet(5.fc))-(l+'"'/^ (14) 

= 9ok + --, — ( hgofe-), (15) 

7**^ is the inverse matrix to gik- We are here particularly interested in the case of dust 
iw = 0, which is an approximation accurate enough to illustrate the method, when 
we integrate the null geodesies back to the surface of last scattering of the CMB 
photons. In the dust case the fluid is not accelerated, and g^i is time- independent. 
goi is not a gradient, when vorticity (the antisymmetric part in Uk) is present. We 
may however use Eqn (11) to transform the symmetric part of Uk to zero. For more 
information we have to look at the field equations. 

5. FIELD EQUATIONS 

The Einstein field equations cannot be integrated as easily as the matter conser- 
vation equations, so we use Taylor expansion again. It is straightforward to write 
them down for the metric (l)-(3) and the matter tensor (13). To zero order one 
obtains, again assuming dust and with po &s matter density on L: 



,fl ,1 ^kl^kl 

a 2. a* 



o I ^kl^kl .,„N 



= -Z-ri 5 ) H T^ihlk — Ikll) -\ J {uJlr{hlkr — hrkl) + UJlki — ^hlrr + hrrl)) 

2 at a^ a^ a* 

(17) 



/d a^ 1 

^a+V-2' 



Skiiz + 2— - o'^Po) 



1 1 

T^kr^lr i Ty^klrr i flrrkl '^rlkr '^krlr) 

a* a* 

+ ^— g (2/lsrr — h„s){hskl + hslk — hkls) 
— J—^ihsrl + hslr — hrls){hskr + hgrk — hkrs)- (18) 

We have written the FRW terms on the Ihs and the terms representing external 
sources on the rhs. We have also kept all nonlinear terms to see the structure of 
the equations, even if these terms can be neglected for a first calculation of the 
effects of external gravitational fields. The equations show that corrections to the 
Friedman equations occur already at zero order, i.e., on the observer world-line L. 
The time dependence of the expansion coefficients is only partly determined by the 
field equations at the given level, thus one has to take higher levels into account. 
Some equations become algebraic relations between the coefficients. We note that 
the restrictions for the expansion coefficients in (l)-(3) found in this way are only 
necessary conditions for g^j^i, to represent a solution of the field equation. Our 
expansion procedure is as yet not based on a well-founded initial value problem, 
and it is therefore not sure that every approximate expressions for g^^i, can be 
extended to a full solution of the Einstein field equations. 

6. THE CMB: ORIGIN OF A COSMOLOGICAL DIPOLE COMPO- 
NENT 

The CMB anisotropy in our local cosmological model can be calculated using a well- 
known relation derived by Panck (1986), Russ et al. (1993) and Dunsby (1997). 
We write the equation as 

where p^ — ^- is the tangential vector of a null geodesic connecting the emission 
(at time ti) of the radiation with the recording event (at time ig); /^(r) the density 
of the radiation field and E = — p^V^ the photon energy relative to V^. Using 
this relation requires integration of nuU-gcodesics, i.e. solution of the geodesic 
equation pf^^ = 0. This can easily be done within the formalism, if one considers the 
additional external terms as perturbations to the known photon motion in a FRW 
geometry. The external or super-horizon modes under discussion affect the lowest- 
order multipoles. As discussed by Bunn et al. (1996), they must be distinguished 
from the familiar statistical fluctuations in T, which correspond to an isotropic 
Gaussian random field and are believed to be generated by quantum processes in 
an inflationary stage. We treat as example the possible presence of an intrinsic 
dipole component {^)dipoie = /""-'^f^ of flic CMB, which is independent of the 
usual kinematic component due to observer motion relative to the CMB rest frame 
(Langlois 1966, Lineweaver 1966). In principle, all terms in (19) may contain dipole 



contributions. The first term results from a spatial gradient of the radiation energy 
density along the past null geodesies and requires the use of perturbed geodesies to 
obtain a nonzero result. The second term is produced by large-scale gravitational 
potential variations and is calculated using the FRW approximation for the null 
geodesies. Its dipole contribution can be written 

with Ai ~ Jf° dt/a. The last term proportional to the acceleration of the fluid is 
absent in the case of noninteracting dust. 

7. CONCLUDING REMARK 

The MAP and PLANCK SURVEYOR satellite missions are expected to measure - 
together with other observations - many cosmological parameters such as the am- 
plitudes and spectral indices of scalar and tensor fluctuations and the densities of 
various mass components of the universe. We emphasize that to these fundamental 
cosmological quantities one has to add further geometrical and kinematical param- 
eters - the spatial curvature of the FRW models is only the simplest one. 
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